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1  INTRODUCTION 


Seismic  waves  travelling  over  regional  distances  (100-1000  km)  from  shal¬ 
low  sources  must  travel  through  the  crustal  waveguide.  There  have  been 
many  studies  of  this  type  of  propagation,  and  a  number  of  attempts  to  ex¬ 
plain  its  behavior  in  terms  of  models  of  earth  structure.  Some  explanations 
have  concentrated  on  the  details  of  horizontally  stratified  wavespeed  struc¬ 
ture,  with  the  wave  equation  being  solved  in  terms  of  normal  modes.1  Other 
explanations  have  used  ray  theory  as  a  starting  point.2  Most  of  these  analy¬ 
ses  have  modelled  earth  structure  in  a  deterministic  manner;  that  is,  models 
were  presented  that  defined  the  seismic  wavespeed  exactly  for  every  point  in 
space. 

For  seismic  waves  above  1  Hz,  a  number  of  recent  studies  have  analysed 
waveform  variability  in  terms  of  a  statistical  model  of  earth  heterogeneities. 
Teleseismic  P-wave  amplitude  and  arrival-time  fluctuations  have  been  ex¬ 
plained  on  the  basis  of  heterogeneities  in  the  lithosphere  and  upper  mantle; 
these  studies  were  sensitive  to  heterogeneities  with  rms  wavespeed  variations 
of  a  few  percent  and  horizontal  scale  lengths  from  a  few  kilometers  to  tens 
of  kilometers.3  These  teleseismic  waves  come  up  toward  the  receiver  stations 
with  angles  not  to  far  from  the  zenith.  (See  Figure  1.)  Propagation  over  re¬ 
gional  distances  from  deep  earthquakes  (well  below  the  crust)  has  also  been 
analysed  in  terms  of  heterogeneities,  and  similar  results  have  been  obtained.4 

In  both  the  above  studies,  the  model  of  deterministic  propagation,  before 
the  addition  of  heterogeneities,  was  a  very  simple  one;  it  consisted  of  incident 
straight-line  rays,  or  plane  waves. 

It  is  the  purpose  of  this  study  to  analyse  regional  propagation  from  shal¬ 
low  sources  from  the  same  point  of  view  as  the  above  studies.  For  this  purpose 
we  must  make  a  simple  model  of  deterministic  propagation  in  the  absence  of 
heterogeneities;  this  is  done  in  Section  2.  Then  the  effects  of  heterogeneities 
are  added  in  Section  3.  We  will  see  that  this  statistical  approach  explains 
many  of  the  features  of  regional  propagation. 

Our  simple  model  consists  of  a  constant-depth  crustal  waveguide  with  a 
constant  wavespeed,  overlying  a  uniform  half-space  with  a  constant  (higher) 
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Figure  1.  Schematic  of  the  arrival  of  seismic  waves  at  a  receiver  on  the  earth's  surface,  for  a  simple 

deterministic  model  of  earth  structure  involving  the  crustal  waveguide.  Waves  from  near  the 
zenith  come  nearly  straight  through;  they  represent  teleseismic  arrivals.  Waves  from  sources 
that  are  below  the  crustal  waveguide  are  refracted  slightly,  but  still  have  a  single,  well-defined 
ray  as  a  deterministic  model  of  the  propagation.  Waves  from  sources  within  the  crustal 
waveguide  can  reflect  many  times  on  their  way  to  a  receiver. 
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wavespeed.  A  real  region  of  the  crust  will  be  more  complicated  than  our 
simple  model,  even  disregarding  small-scale  heterogeneities.  For  example, 
the  depth  of  the  waveguide  may  be  variable  with  geographical  position,  or  the 
wavespeed  may  have  a  small  gradient,  or  a  series  of  step  changes  within  the 
waveguide.  We  will  see  that  when  we  add  the  heterogeneities,  the  complexity 
introduced  into  the  received  waveforms  is  so  drastic,  that  if  we  had  started 
with  a  more  complicated  deterministic  model,  it  would  not  have  qualitatively 
changed  our  conclusions. 

In  this  study  we  will  ignore  the  problem  of  conversion  between  compres- 
sional  and  shear  waves.  Our  study  may  be  thought  of  as  initially  applying  to 
SH  propagation,  since  there  is  no  conversion  of  SH  waves  in  our  deterministic 
model.  Again,  the  variations  introduced  by  the  heterogeneities  are  so  signif¬ 
icant,  even  for  pure  SH  waves,  that  our  qualitative  conclusions  will  probably 
apply  in  general,  but  this  is  clearly  an  area  deserving  of  further  first-order 
study. 
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2  CRUSTAL  WAVEGUIDE  PROPAGATION 


Deterministic  propagation  in  the  crustal  waveguide  at  high  frequencies 
can  be  analysed  from  the  point  of  view  of  ray  theory.  Our  model  will  be  a 
uniform-speed  waveguide  of  thickness  h  and  wavespeed  /?,  with  a  uniform- 
speed  half-space  underlying  it.  The  speed  in  the  underlying  medium  will 
be  designated  /?2.  We  will  take  a  source  at  depth  d  at  a  distance  R  from  a 
receiver  station  on  the  surface.  (See  Figure  2.) 

The  solution  to  our  simple  model  is  straightforward,  and  consists  of  a 
sequence  of  straight-line  rays  that  are  reflected  n  times  from  the  surface 
or  bottom  of  the  waveguide,  where  n  is  0,  ±1,  ±2  ...up  to  a  maximum 
determined  by  critical  reflection  at  the  bottom  of  the  waveguide.  This  critical 
reflection  angle  is  determined  by  the  ratio  of  the  wavespeeds  in  the  two  regions 
of  the  earth: 

cos0critical  <  W2. 

Let  ©„  be  the  angle  of  the  nth  ray  with  the  horizontal.  The  requirement  that 
0n  <  €)critical  can  be  more  usefully  expressed  by  the  equation: 

\d  +  2nh\  <  [(&//?)2-l }1/2R. 

The  number  of  rays  N  can  be  expressed  approximately  from  this  equation: 
N  =  {R/h)[W/3f-l]''2. 


The  intensity  of  each  ray  can  be  calculated  from  simple  spherical  spread¬ 
ing,  because  each  ray  can  be  unwrapped  from  its  reflections  into  a  simple 
straight  line  having  a  length  of  R/cosQn.  Therefore,  the  intensity  of  the  nth 
ray  is: 

In  =  S(Qn)cos2ejR 2  =  S(Qn)/[R2  +  (d  +  2  nh)2) 

where  S(0„)  is  the  angular  distribution  at  the  source.  From  here  on  we  will 
assume  for  simplicity  that  the  source  is  isotropic;  that  is,  ^(©n)  =  1.  If 
we  normalize  the  intensity  pattern  by  having  the  first  ray  (with  n=0)  have 
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Figure  2.  Diagram  of  rays  from  a  source  in  the  crustal  waveguide  to  a  receiver  on  the  surface.  Each  ray 
has  an  angle  8n  with  the  horizontal.  The  wavespeeds  in  the  crustal  waveguide  and  below  are 
0  and  02.  respectively.  Rays  with  angles  above  the  critical  angle  get  lost  in  the  deep  earth. 


an  intensity  of  unity,  then  it  is  easy  to  show  that  the  last  ray  will  have  the 
minimum  intensity,  and  it  will  be 

'min  =  (W*)2 

The  arrival  time  of  each  ray  is  easily  calculated  from  the  distance  travelled 
( R/cosQn ): 

K  =(«//3){l  +  ((<*  +  2nh)IR?}''\ 

Note  that  if  d  <<  h,  then  the  two  rays  corresponding  to  positive  and  negative 
n  come  in  very  close  together;  there  is  a  twinning  of  arrivals  associated  with 
the  source  and  its  image  reflected  in  the  earth’s  surface.  The  patterns  that 
result  for  several  ranges  of  propagation  in  a  crustal  waveguide  with  realistic 
average  values  for  the  various  characteristic  parameters  are  shown  in  Figure 
3. 

It  is  useful  to  calculate  the  approximate  spacing  between  the  discrete 
arrivals  for  later  comparison  with  the  effects  of  heterogeneities.  Let  A tn  be 
the  spacing  between  the  nth  and  (n+l)th  rays: 

Af„  =  (2ft//))f„/[l  +  (T2 

where 

£n  =  \(d  +  2nh)/R\. 

This  expression  does  not  account  for  the  twinning  due  to  positive  and  neg¬ 
ative  n.  It  is  of  interest  to  note  that  near  the  first  arrival,  the  separation  is 
much  smaller  than  near  the  last  arrival: 

At,  =  {2h/P)(2h/R) 

Atmax  =  (2  h//3). 

For  example,  if  h=30  km,  R=600  km,  and  /?  =3  km/s,  then  At,  =  2  s 
and  Atmax  —  20  s. 
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UNCLASSIFIED 

Figure  3.  (U)  Examples  of  expected  high-frequency  waveforms  for  our  simple 

deterministic  crustal  waveguide  parameters,  with  no  heterogeneities. 

The  plots  are  for  two  different  ranges  (300  km  and  1000  km)  and  two 
different  source  depths  (2  km  and  10  km).  The  parameters  of  the  crust  are: 
h=30  km,  *3=3  km/s,  and  02=5  km/s. 
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3  THE  EFFECTS  OF  HETEROGENEITIES 


The  simple  arrival  structure  described  in  Section  2  is  not  observed  exper¬ 
imentally.  The  arrival  pattern  is  more  complex,  with  intensity  as  a  function 
of  time  being  an  apparently  random  jumble  that  begins  about  where  the  first 
arrival  is  expected,  and  continues  until  dying  away  in  the  noise  many  tens  of 
seconds  later. 

We  will  calculate  the  effect  of  small-scale  heterogeneities  on  our  determin¬ 
istic  arrival  structure  by  use  of  the  path-integral  theory  of  wave  propagation 
through  random  media,  developed  in  large  part  for  the  case  of  acoustic  prop¬ 
agation  in  the  ocean.5-11  To  begin  with,  we  need  a  statistical  model  for  the 
heterogeneities.  For  the  purpose  of  this  illustrative  study,  we  will  take  an 
anisotropic  Gaussian  correlation  function  with  scale  lengths  in  the  horizon¬ 
tal  (Lv')  and  vertical  (L//)  and  an  rms  fractional  variation  in  wavespeed  of 

/i. 


We  emphasize  that  we  are  taking  a  simplistic  picture  of  the  medium  in  as¬ 
suming  that  it  is  composed  of  the  superposition  of  two  disparate  structures: 
First,  a  very  simple  deterministic  structure  (a  plane  layer  over  a  homoge¬ 
neous  half-space),  and  second,  small-scale  random  velocity  perturbations  in 
the  crustal  layer.  Other  effects,  which  are  certainly  of  importance  in  specific 
regions,  might  be  topography  on  the  surfaces  of  the  crustal  layer,  or  reflect¬ 
ing  surfaces  within  the  crust.  The  point  of  our  study  is  that  much  of  what 
is  actually  observed  may  often  be  explained  by  random  volume  velocity  per¬ 
turbations,  so  that  such  perturbations  should  not  be  ignored  in  a  first-order 
explanation  of  observations,  particularly  at  high  frequency. 

We  will  take  the  values  of  velocity-fluctuation  quantities  from  the  stud¬ 
ies  that  have  been  carried  out  on  teleseismic  waves  arriving  at  large  arrays, 
and  from  recent  studies  of  regional  propagation  from  deep  earthquakes.  The 
relation  between  crustal- waveguide  propagation  and  teleseismic  effects  is  de¬ 
pendent  on  the  anisotropy  of  the  medium  fluctuations,  which  is  not  well 
known.  We  are  only  giving  order-of-magnitude  examples  here;  in  the  study 
of  a  specific  geographical  region,  careful  consideration  to  specific  anisotropy 
models  should  be  given. 
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Typical  values  of  these  fluctuation  quantities,  from  the  teleseismic  study 
and  the  study  of  regional  propagation  from  deep  earthquakes,  might  be: 
fx  ~  2%,  Lv  ~  5  km,  and  Lh  ~  10  km.  However,  the  results  are  quite 
sensitive  to  the  exact  choices  of  values,  and  to  the  particular  spectrum  chosen. 
Therefore,  a  more  quantitative  study  must  be  done  to  compare  with  actual 
data  from  a  specific  geographical  region. 

3.1  Fluctuation  Regime 


Various  regimes  of  wave  fluctuation  behavior  were  identified  in  Reference 
5  and  6,  and  quantities  were  defined  that  are  necessary  for  determining  the 
regime  in  which  a  particular  situation  lies.  These  quantities  are: 


•  The  strength  parame*~r  $,  which  is  the  rms  phase  that  a  wave  would 
experience  if  it  travelled  exactly  along  the  deterministic  ray.  This  would 
be  the  actual  rms  phase  only  if  the  heterogeneities  were  weak  enough, 
but  even  when  the  heterogeneities  are  strong,  the  calculation  of  the 
strength  parameter  is  crucial  to  the  wave-fluctuation  analysis. 

•  The  diffraction  parameter  A,  which  is  the  square  of  the  ratio  between 
the  size  of  the  Fresnel  zone  and  the  transverse  correlation  length  of  the 
medium,  averaged  over  the  deterministic  ray.  Small  diffraction  param¬ 
eter  means  that  diffraction  is  small,  so  that  the  concepts  of  generalized 
geometrical  optics  are  quite  useful.  Large  diffraction  parameter  means 
that  diffraction  is  very  important;  in  this  situation  concepts  of  uncor¬ 
related  multiple  scatter  are  probably  useful. 

•  The  various  regimes  in  A  —  $  space  are  shown  in  Figure  4.  We  need 
to  define  a  few  more  quantities  in  order  to  calculate  these  parameters. 
They  are: 

•  The  wavenumber  of  the  propagating  seismic  wave,  which  we  will  denote 
by  q.  For  example,  a  1-Hz  shear  wave  travelling  in  a  crust  with  a  speed 
of  3.14  km/s  has  a  wavenumber  q  =  2  km-1. 
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Figure  4.  A-4>  diagram  showing  the  different  regimes  of  wave  fluctuation  behavior  for  waves 
propagating  through  random  media  (from  Reference  6). 


11 


•  The  parallel  (Lp)  and  transverse  (Lx)  correlation  length  of  the  medium, 
where  parallel  means  along  the  deterministic  ray,  and  transverse  means 
perpendicular  to  the  deterministic  ray.  In  an  anisotropic  medium,  these 
can  often  be  expressed  in  terms  of  Ly  and  Lj/. 


In  terms  of  defined  quantities,  the  following  results  for  the  strength  and 
diffraction  parameters  follow: 

$2  =  q2ti2LPR  =  16 

A  =  R/(GqL2T)  =  3.3 

A$2  =  ( l/6)qn2R2LP/L2T  =  53 

A$  =  (1  /6)fxR(LPR)1/2/Ll  =  13 

and  putting  in  our  so-called  typical  numbers  yields  the  numerical  values 
indicated,  for  a  propagation  range  of  1000  km. 

The  fact  that  both  A$  and  $  are  greater  than  unity  tells  us  that  we  are 
in  full  saturation.  In  this  regime  we  have  many  uncorrelated  microrays;  that 
is,  each  deterministic  ray  is  broken  up  by  the  heterogeneities  into  many  rays. 
This  remains  true  as  one  moves  to  higher  frequency. 


3.2  Intensity  Fluctuations 


In  the  fully  saturated  regime,  the  intensity  of  a  given  ray,  at  least  at 
a  single  frequency,  is  Rayleigh  distributed.  That  means  that  the  intensity 
has  a  simple  exponential  distribution  function.  As  long  as  the  analysis  is 
being  done  on  a  limited  frequency  band  that  does  not  exceed  the  coherent 
bandwidth,  the  result  is  that  each  of  the  arrivals  in,  for  example,  Figure  3, 
will  have  an  intensity  taken  from  this  exponential  distribution  function,  with 
mean  value  given  by  the  deterministic  result.  However,  as  will  be  shown 
below,  the  case  of  interest  here  has  a  very  small  coherent  bandwidth.  In 
that  case,  the  observation  can  be  thought  of  as  being  within  a  small  interval 
of  time,  so  small  that  the  time  spread  of  the  received  pulse  is  much  larger 
than  the  observation  interval,  and  the  intensity  within  each  observation  time 
interval  will  be  a  random  variable  with  a  large  variance. 
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3.3  Time  Spread 


In  the  fully  saturated  regime,  a  delta-function  pulse  emitted  by  the  source 
is  spread  in  time  by  an  amount  given  by5 

t.  =  a*V(«J){i  +  £}. 

We  can  compare  this  spread  to  the  separation  between  the  deterministic 
arrivals  by  taking  the  ratio  of  the  two: 

r„/A(n  =  A*V(2?A)!l+a3/2/«n- 

The  £n  at  small  n  are  small  compared  with  unity  (on  the  order  of  h/R), 
and  rise  to  a  maximum  of 

£max  =  [(/W/?)2  “  l] 

A  few  simple  substitutions  of  numerical  values  reveals  that  the  ratio  of 
t/ At  is  always  greater  than  unity.  Thus  the  spread  due  to  heterogeneities  is 
always  larger  than  the  separation  between  deterministic  arrivals,  resulting  in 
a  waveform  that  essentially  appears  to  be  random. 

For  example,  for  our  typical  numbers  at  1000  km  range,  we  find  that 

Ti  =  9  s 
tisi  —  25s 

and  we  see  that  in  all  cases  the  spread  fills  in  any  gaps  one  might  have 
observed  in  the  deterministic  propagation. 

The  coherent  bandwidth  of  the  propagation  is  the  inverse  of  the  time 
spread.  The  numbers  above  imply  that  the  coherent  bandwidth  is  in  the 
range  of  0.04-0.10  Hz. 
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3.4  Angular  Spread 


In  the  fully  saturated  regime,  the  energy  that  starts  from  a  point  source 
and  arrives  at  a  point  receiver  has  travelled  within  a  spatial  region  that  is 
determined  by  the  medium  fluctuations  rather  than  the  finite  wavelength  of 
the  seismic  wave.  The  size  of  this  spatial  region  also  determines  the  spread 
in  arrival  angle  of  the  energy  as  it  arrives  at  the  receiver. 

The  determination  of  this  angular  spread  is  dependent  on  a  function 
known  as  the  phase  structure  function,  which  is  the  variance  of  phase  dif¬ 
ference  between  two  nearby  rays  that  start  from  the  source  and  end  near 
the  receiver  with  a  small  separation  x.  In  the  case  of  a  Gaussian  correlation 
function,  the  phase  structure  function  D  is  given  by 

D(x)  =  V(zlLT)\ 


The  rms  angular  spread  is  determined  by  the  value  of  x  (call  it  x0)  at 
which  Z)(x0)=l.  Thus: 

x o  = 

and 

Grms  =  l/(xaq)  =  $/(qLt). 

For  our  typical  case,  in  which  $  =  4,  <7  =  2  km-1,  and  Lj  =  5  km,  we 
have  Orms  =  0.4  radians,  or  23°.  Note  that  the  angular  spread  scales  with 
the  square  root  of  range,  and  is  independent  of  frequency. 

Thus,  an  array  that  attempts  to  determine  the  location  of  the  source  by 
beamforming  will  find  that  the  inherent  inaccuracy  due  to  the  heterogeneities 
has  a  standard  deviation  on  the  order  of  20°  at  1000  km. 
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3.5  Averaging  and  Calibration  Distances 


The  distance  xQ  represents  the  distance  over  which  the  arriving  wave 
field  is  decorrelated.  Therefore,  if  one  wants  to  go  back  and  verify  the  exact 
waveform  received  in  a  particular  event  by  repeating  the  transmission  by 
some  means,  then  one  must  place  both  the  source  and  the  receiver  within  a 
fraction  of  x0  of  the  original  transmission.  This  would  be  true  for  both  the 
horizontal  and  vertical  position.  In  our  typical  case,  the  value  of  xa  is  about 
1  km. 

By  the  same  token,  if  one  were  attempting  to  average  out  the  effects 
of  heterogeneities,  one  would  get  approximately  independent  samples  every 
x0i  and  therefore  many  receivers  spaced  by  that  distance  would  give  a  good 
average. 

However,  this  conclusion  is  really  not  accurate  for  media  that  are  mul¬ 
tiscale  rather  than  Gaussian.  We  would  expect  that  the  larger  averaging 
region  would  be  affected  by  larger-scale  medium  variations,  and  the  variance 
would  not  go  down  proportionally  to  the  number  of  receiving  stations.  This 
kind  of  calculation  can  be  done  for  any  given  model  of  medium  fluctuations. 
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4  DISCUSSION  AND  CONCLUSIONS 


An  intuitive  understanding  of  regional  propagation  of  high  frequency  seis¬ 
mic  waves  for  sources  that  lie  within  the  crustal  waveguide  has  been  presented 
in  terms  of  a  very  simple  deterministic  propagation  model  combined  with  the 
effects  of  small-scale  heterogeneities.  Our  simple  deterministic  model  implies 
that  the  underlying  pattern  of  arrivals  is  a  series  of  sharp  arrivals  associated 
with  the  multiray  structure  expected  in  a  waveguide.  Reasonable  numer¬ 
ical  values  for  heterogeneities  in  the  crust  yield  very  significant  effects  on 
the  received  arrivals.  These  effects  include  variations  in  amplitude,  arrival 
time,  and  arrival  angle.  As  a  result  of  these  fluctuations,  the  received  wave¬ 
forms  will  lose  their  detailed  connection  with  the  underlying  arrival  pattern, 
and  will  become  a  jumble  of  arrivals  whose  spread  in  time  overlaps  with 
the  nearby  deterministic  arrivals.  The  variations  in  intensity  of  each  arrival 
caused  by  the  heterogeneities  are  large. 

The  angular  spread  of  the  arriving  energy  can  be  calculated,  and  for  1000- 
km  propagation  has  an  rms  value  of  20°  or  more.  This  angular  spread  also 
implies  that  nearby  receivers  will  not  be  coherent  with  each  other  at  distances 
comparable  to  1  km,  and  that  the  use  of  calibration  shots  of  some  kind  in 
order  to  measure  the  propagation  from  a  particular  source  to  a  particular 
receiver  must  be  done  under  the  limitation  that  the  calibration  shot  must  be 
within  that  1-km  distance. 
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